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We perform a global fit to Higgs signal-strength data in the context of light stops in Nat- 
ural SUSY. In this case, the Wilson coefficients of the higher dimensional operators medi- 
ating g g ^ h and /i — ^ 7 7, given by Cg, c^, are related by Cg = 3 (1 + 3 ag/ (2 7r))c^/8. 
We examine this predictive scenario in detail, combining Higgs signal-strength constraints 
with recent precision measurements of mw, Br{B — )• Xg 7) constraints and direct collider 
bounds on weak scale SUSY, finding regions of parameter space that aie consistent with all 
of these constraints. However it is challenging for the allowed parameter space to reproduce 
the observed Higgs mass value with sub-TeV stops. We discuss some of the direct stop 
discovery prospects and show how Higgs search data can be used to exclude light stop pa- 
rameter space difficult to probe by direct collider searches. We determine the current status 
of such indirect exclusions and estimate their reach by the end of the 8 TeV LHC run. 



I. INTRODUCTION 

The discovery of the Standard Model (SM) Higgs (assuming that the discovered boson is the 
SM Higgs) puts Supersymmetry (SUSY) to a new test. The possibility that SUSY solves the hierar- 
chy problem without re-introducing fine-tuning i.e. the paradigm of Natural SUSY (NSUSY) [[T]- 
[3 can now be examined in the light of increasingly precise measurements of the properties of the 
Higgs. 

NSUSY predicts new particles near the electroweak (EW) scale, which necessarily must affect 
the properties of the Higgs if they are to stabilize this scale. The minimalistic scenario for NSUSY 
focuses on the vestige of the SUSY spectrum which is required to be light in order to keep the 
fine-tuning of the theory reasonably small. In this case, the most significant impact of new states 
on Higgs phenomenology is through the presence of light stop states. This scenario is not just 



2 



motivated by simplicity, but also by the lack of evidence for SUSY to date, indicating that a weak 
scale SUSY spectrum needs to be non-generic to satisfy collider constraints. The states directly 
related to naturalness (primarily the stop and Higgsinos) are especially challenging, and model 
independent collider bounds are weak or non-existent. 

Conversely, the study of the impact of an NSUSY scenario on the properties of the Higgs 
benefits from the enormous effort expended by the experimental collaborations in refining the 
accuracy and precision of the reported Higgs signal strength measurements. When considering the 
experimentally resolvable impact of NSUSY, indirect probes, e.g. through a fit to Higgs properties, 
electroweak precision data and flavor physics may well be more powerful in constraining many 
minimal scenarios than direct searches for some time. This is the line of reasoning we develop 
in this paper, where we examine the current constraints on minimal NSUSY from these indirect 
probes. 

The outline of this paper is as follows. We briefly review and introduce NSUSY in Section |n} 



In Section III we then review the impact of the NSUSY spectra on the properties of the SM Higgs 



through modifications in the loop-induced h g g and /i — 7 7 couplings. Further, in Section IV 
we work through the constraints set on NSUSY from global fits to Higgs properties in the presence 
of light stops. We advance such studies by using a more complete global fit (now including 48 sig- 
nal strength channels, including ICHEP and post ICHEP data updates)^. We make use of the fixed 
relationship between the Wilson coefficients (including the QCD matching corrections) for /i77 
and hgg in the case of NSUSY to perform a one -parameter fit and then directly map the allowed 
parameter space in global Higgs fits into the allowed stop space. Further, we determine 95% 
confidence level (C.L.) exclusion limits on the stop parameter space derived from Higgs search 
data. We then consider constraints from BR(i? — Xs 7) and recent precision measurements of 
mw at the Tevatron, as these results, which are under excellent theoretical control, are sensitive 
to light stops. The (statistically insignificant, but interesting) deviations from the SM predictions 
in these observables in a weak scale NSUSY scenario could offer some further resolution on the 
allowed stop parameter space, if NSUSY exists. Ascribing these deviations to the effect of stops in 
NSUSY, such stop states are consistent with the results of the global fits to Higgs signal strengths, 
as we will show. Finally, we also take into account direct collider bounds. In SectionlVlwe discuss 



^ Some past papers that have consistently examined earlier versions of the Higgs dataset in this context, with varying 
degrees of sophistication, are lIs mTl . 



3 



Field 


Spin 


SU(3)c X SU(2)l X U(1)y 


Ql = (tLM) 





(3,2,1/6) 


f* 

iX 





(3,1,-2/3) 







(1,2, +1/2) 


Hd={HlH-) 





(1,2,-1/2) 




1/2 


(1,2, +1/2) 




1/2 


(1,2,-1/2) 


9 


1/2 


(8,1,0) 



TABLE I: The minimal NSUSY field content. In our analysis, gluinos and a heavy linear combination of the 
Higgs doublets will be further integrated out. 

the interplay of these constraints and determine the allowed parameter space that remains. We in- 
clude the limits coming from the Higgs mass measurement and estimate the degree of fine-tuning 
incurred, finding that some regions of the stop parameter space that give a reasonable global fit 



also allow the Higgs mass to be raised to its observed value. In Section VI we discuss the cur- 
rent exclusion bounds that can be derived using these indirect probes of stop parameter space. To 
study the future prospects of such limits by the end of 2012, assuming that the experimental error 
in Higgs signal- strength measurements scales down as ~ 1/y/Cint, we consider two hypothetical 
cases: 1) the current pattern of best-fit signal-strength values does not change, and 2) the dataset 



evolves to converge on the SM expected signal strengths. Finally, in Section VII we conclude. 



II. NATURAL SUSY 



Naively, in generic SUSY scenarios motivated as a solution to the hierarchy problem, one 
expects all superpartners near the electroweak scale, with the soft breaking mass scale Msusy not 
higher than 0{1) TeV The experimental picture emerging from the LHC is in growing tension 
with this expectation. After searching in many typical discovery channels, and reaching a peak 
sensitivity of roughly C(1.5TeV)/ C(100)'s GeV, for coloured/electroweak SUSY states Ml, no 
statistically significant experimental excess has been reported to date. On the other hand, to avoid 
destabilizing the electroweak scale when Msusy ^ v (without fine-tuning), only a minimal set of 
SUSY particles have to be light (< 1 — 2 TeV) [[I]-|7l. The stop soft masses are directly connected 
(at one-loop) to the Higgs mass scale (or Z mass) through the sizable top coupling, so fine-tuning 
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considerations require them to be light. Although sbottoms (6) do not directly affect the fine-tuning 
of the Z mass, 6^ is required to be light as it is linked to the tL mass scale by SU(2)l symmetry. 
The Higgsino mass coming from the /i term^ is tied to tree-level contributions to the Higgs mass. 
With gaugino masses assumed heavy, of order Msusy, light charginos and neutralinos are almost 
pure Higgsinos, with mass given by fi up to corrections of order v/Msusy that we neglect. To 
a lesser degree, the gluinos g are also expected to be light due to their contribution to two-loop 
corrections to the Higgs mass parameter, which leads to the rough estimate nig ^ 2m^ (4mj) for a 
Majorana (Dirac) gluino [6J. In practice, we can also decouple these somewhat heavier gluinos in 
our analysis. 

The states just discussed are listed in Table |I} These are the states whose impact on Higgs 
signal strengths and low-energy precision measurements we will focus on in this paper. We will 
not consider light staus (f^*"), which can also modify the Higgs decays to photons, see Refs. [|8l|9]l 
for recent studies. We neglect these states as we assume a moderate value of tan (3 ^ 0(10), for 
which effects are negligible compared to the stop contribution. We also neglect the effects of 
the b on the Higgs mass and in the loop corrections to the Higgs signal strength parameters for the 
same reason. 

Regarding the SUSY Higgs sector, we will consider the decoupling regime in which only one 
Higgs doublet remains hght^, while the second doublet, with mass controlled by the pseudoscalar 
mass iriA, has a mass ~ TeV. In this limit, the couplings of the lightest Higgs h to fermions and 
gauge bosons approach their SM values, and we will only consider deviations in the loop-induced 
couplings of the light h to photons and gluons (see next section). 

In order to fix our notation, we write now the parts of the low-energy Lagrangian most rele- 
vant for our analysis. This Lagrangian is not supersymmetric as it applies below the scale of the 
heavy SUSY particles (with masses ~ Msusv ~ 1 TeV). Supersymmetric relations between some 
couplings are broken and one should introduce different couplings, to be matched to the supersym- 
metric theory at the scale Msusv- In practice, the hierarchy between Msusv and the electroweak 
scale is mild and one can neglect most of these breaking effects. One important exception is the 

^ In this paper we distinguish signal-strengths with a subscript, /i^ for a final state i, from the /i parameter in NSUSY, 
which carries no subscript. 

^ This choice is supported by the results of Ref. [171, where the 2HDM is studied in the light of the Post-Moriond 
Higgs data and no compelUng region in the parameter space of the MSSM consistent with the data was identified, 
except the decoupling limit. 
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Higgs quartic coupling, as will be discussed later on. 

From the superpotential W = ^H^ ■ Hu + htQ l ■ H2t% the Lagrangian gets the terms 

5C = fiHd ■ Hu + KQl ■ Hufk + htQL ■ HutR + ^Ql ■ Hui*R + h.c. 



(1) 



where ■ stands for the SU(2)l product: H^ ■ Hu = H^Hu — H^H^^^ and we have suppressed 
SU(3) indices for (s)quarks. (Weyl) fermionic fields are contracted in the usual way: H^H^ = 
i{H^Y a2Hu- Eq. ([T]) contains a Dirac mass term /i for Higgsinos, the top Yukawa coupling 
between top quarks and the Higgs, and the related Higgsino-squark-quark couplings, relevant for 
the contribution of Higgsino-stops to Br (5 — )■ Xs'y). 

The scalar potential for Higgses and squarks is well known and includes supersymmetric F 
and D terms and soft-SUSY-breaking terms. We can always perform a rotation of the full Higgs 
doublets Hu d to the doublets Hi h. 




(2) 



cos /3 sin (3 
— sin (3 cos (3 

where tan (3 = / {Hf) = Vu/vd and H^ denotes the SU(2) conjugate, H^ = -ia'^H^ = 
(—H^,H^*), in such a way that Hi is the doublet involved in electroweak symmetry breaking 
while Hh does not take a vacuum expectation value. In the Higgs decoupling limit, with large 
pseudoscalar mass, Hh is in fact composed of the heavy fields H^, A^, H^, while Hi is SM-like 
and contains the light Higgs h and the Goldstones. The quartic Hi coupling determines the light 
Higgs mass as usual and is the prime example of a coupling that receives sizeable SUSY-breaking 
corrections (that help in increasing the Higgs mass above its tree level minimal SUSY value below 
mz). Such corrections will be discussed in Subsection|VjA. 
Finally, the stop masses are given by the mass matrix 



Mi 



'mi. 


Mln 






) C2/3 mt(At + /i/tan/3) 


Ml, 


Mln_ 




mt{At + yu/ tan/3) 





(3) 



where we used C2p = cos 2/3, rrit is the top mass and Mq^, M^^, At are soft SUSY-breaking 
masses. The stop mixing angle 9t relates the interaction eigenstates iL,R to the mass eigenstates 
1 1,2 by the rotation 



cos 6t — sin 6f 
sin Of: cos 6^ 



(4) 
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The mixing angle Of is taken in the interval (— 7r/2, 7r/2) and defined by 

L/JX 

with the signs of M^j^ — and Mlj^ determining the quadrant of 29^. With this definition of 
9^, one automatically guarantees m^^ < m^^, with 



cos 29i = ^' , sin 29i = , (5) 



ml - ml ^ {5mf = ^ {Ml, - Ml^f + 4M|^ . (6) 

Finally, neglecting sbottom mixing (proportional to rrib), the light sbottom has mass m? = M? + 
ml — m|(l/2 — cos 2/3, while the heavy hn is decoupled. 



III. LOOP-LEVEL CORRECTIONS TO HIGGS PROPERTIES IN NSUSY 

In this section, we review the loop-level NSUSY corrections to the couplings of h to photons 
and gluons."^ The leading correction from stops to the gluon-fusion process is given by |fT9ll20ll 

(^{99 h) T{h gg) 2 

|1 + rg| , (7) 



a^^igg^h) T^^^ih^gg) 
where 

^ _ Cgias)Fg{mi^,mi^,9i) 
9 F^^^{mt,mb- ■ ■) 

Here Cg{as) is a factor that takes into account higher order QCD corrections- see the discussion 
below- and the Fg functions are defined as follows 

mi 



Fg{mi^,mi^,9i) = ^ . . ^ Fo(r,), (9) 

- - 
i=ti,t2--- 

F/^(m„m,---) = Yl (1 + ^) ^ -2/ (1.41 + 0.14 (10) 

i=t,b- ^ ^ ^ 

where n = ml/{Am^), Fo(r) = [r - /(r)] /r^ and Fi/2(r) = [r + (r - l)/(r)]/r2 with /(r) = 
arcsin^ ^/T for r < 1 while, for r > 1, 



fir) = ~\ 



\ 1 + Vl - . 

log , — ITX 

1- Vl-r-1 



2 



(11) 



^ The effects of NSUSY loop corrections on the decay T{h — > Z7) are neglected as the leading effects in these 
scenarios come from charged scalars in the loop, and we are considering the decoupling limit in the scalar sector 
(consistent with the minimal version of NSUSY). 
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See the Appendix for the SM inputs used in determining the numerical value of F^^ above. 
The couplings g^^i. are given, in the decoupling limit ll2T1l . by 

dHuu"^^ = -1 - UiD^^UU) - ^ [(3-44)Kt,|t^)r + 44|(t,Mr] , (12) 

where = mt{At + jj,/ tan/3) = (m?_^ — m?^) sin cos has been defined in Eq. (3j); the 
{ii\iL,R) can be directly read from Eq. (|4]); and we have included the D-term contributions pro- 
portional to cos 2/3, although their effect is negligible. The sign of Ylii=i 2 9huu ^« negative if 
the the stop sector is dominated by a light eigenstate of pure chirality, and positive if the term 
in dominates, M^^ > 4:mf{m'^^ + m^J. In the no-mixing case, we expect a reduction of 
a{gg —^h).ln the maximal-mixing case, the enhancement depends on the separation between the 
two eigenstates. 

The decay width T{h 77) is also modified by stop loops through the same function in 



Eq. ( 10), as the non-Abelian nature of QCD is irrelevant for the leading-order loop function. One 



finds [19] the correction 



r(/^^77) , |2 _ iV.gf C^(a,)F,(m,-^,m,-„g,-) 



rs^(/i^7 7)~' ' ^ FfA'^(mt,iy,mb 

The SM contribution is given by 

F^^\muW,m^---) = F^{Tw)+ iV, i7^i/2(r,) (l - ^) ^ 1/(0.155 - 0.002 z) , (14) 

i=t,h- 

where Fi{t) = [2 + 3r + 3(2r — l)/(r)] /r^, Nc is the number of colours, and Qi is the elec- 
tric charge with e factored out. The matching correction in this case is given by Cj{as), to be 
further discussed in the next section. 

In the no-mixing limit, the effect of the stop is to increase the branching ratio to photons. 



When the mixing term in Eq.(12) dominates, the effect is to reduce it. The total cross section 
(7(7 — )■ /i — > 77 accounts for those two opposite effects, but the combination tends to follow 
whatever is the sign of the gluon fusion modification. 

IV. CURRENT CONSTRAINTS ON NSUSY 
A. Effective Theory Approach 



It is useful to consider the approximation that all of the light NSUSY states are still heavy 
enough to be integrated out giving local operators. We can then fit to the data directly using the 
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effective Lagrangian^ 



HD 



CWB9192 H^^ajjj^ ^a,u 

2 A2 ^ 



(15) 



where (71, (72, Q'i are the weak hypercharge, SU(2) gauge and SU(3) gauge couplings and the scale 
A corresponds to the mass of the NSUSY states integrated out. The effects in NSUSY appear at 
the loop level, so we find it convenient to rescale the Wilson coefficients as Cj = Cj/{l6n'^). In 
this case, using the results of Ref. ||22]| . the effect of the operators in ([15]) is 



SM 

gg~*h 



,2 7, 



pSM ^2 



) J- ft— >7 7 



ft— >'7 7 



1 v'^ 



pSM A2 



(16) 



Here = cw + cb — cwb- We can translate the effect of stops in the language of local operators 



by inspecting our expressions in Sec. Ill and those in Eq. ( |T6| ), 



til 



GJa. 



A2 



N,QiC^{a,)F, 



g ■ 



(17) 



These relationships are only approximate in the sense that the limit ^ ^ '^z* should be taken 
in the loop functions r^, to match onto the local operators. If the two stop mass eigenstates can 
be integrated out simultaneously, the matching can be directly performed by expanding the loop 
functions in this limit, and one obtains^ 



1 
3 



2 
ti 



m- 4 m- m- 

12 



(18) 



While, in this limit, the QCD matching corrections take the simple form 

25 as ^ / \ Sa^ 



Gg{a,) = 1 + 



6 TT 



1 + 



(19) 



These perturbative corrections are the matching corrections due to top squarks in the loops that do 
not cancel when a ratio is taken with the SM contribution to these loops. This correction factor 



^ We do not include CP violating operators, assuming that all non-SM CP violating phases of the states integrated 
out are negligible. It has been argued that this can be naturally accomplished in NSUSY when additional 
assumptions are employed concerning R-parity for example. Our assumption is conservative as the existence of 
large CP violating phases would only increase the constraints on an NSUSY spectrum. 

^ Here we have neglected D term contributions, although we retain the effect of D terms in some of the numerical 
results presented. These corrections are negligible except for 771^^ < 150 GcV masses. They introduce a (minor) 
tan/3 dependence into the definition of Fg in the local operator approximation when retained. 



9 



is obtained in Ref. [|23l in the limit where gluino effects and the effect of squark mixing was 
neglected.' This approach also neglects running that would sum large logs if a two stage matching 
was employed, integrating out each stop eigenstate in sequence.^ 

Whether one integrates out the stops and matches onto the local operator approximation or 
not, there is a relationship between the NP effects on Cgg^h and (Jh^-^-j that is independent of 
the stop mass parameters in the minimal NSUSY limit. In the local operator approximation, the 
relationship is simply 



where we see how the ratio ~ 3/8 is determined by the stop quantum numbers. This is a con- 
sequence of assuming that the only BSM contribution to both the 77 and gg loops comes from 
stops. This strong relationship will be relaxed in less minimal scenarios. For example, light x^'s 
with mass m^^ ~ fi^ (in the decoupling limit) would in principle also contribute to the 7 7 loops. 
However, the Higgs couples to the higgsino as h and a large mixing between wino and 

higgsino eigenstates would be required. As we are considering the large gaugino mass limit in 
NSUSY, M2 ^ /i, V, this mixing scales as ~ sin^ /3/(M|) and is suppressed, so that we can 
neglect the chargino contribution to h'yy. 

We also utilize this effective Lagrangian to examine the issue of efficiency corrections to the /ij 
when high-dimension operators are present. We find that such efficiency corrections to event rates 
are very small and neglect them. See the Appendix for details. 



In this section we describe our method and results for globally fitting to Higgs signal strength 
data in the scenario discussed above. Here we only briefly review the fit procedure, the details of 
our fit method are given in [|26] - l28ll .^ Our fit incorporates the recently released 7 and 8 TeV LHC 

^ It has been pointed out that rrig 00 leads to mixed stop-gluino UV divergences 1241 requiring extra counter- 
terms, but this technical requirement is not a barrier to the numerical investigations we perform. The full matching 
correction is given in Ref. |25 1: the gluino contributions and stop mixing effects are a small correction to the ~ 5% 
matching correction we consider. 

* There are also perturbative corrections to the matrix element of the local effective operator h Gj^^G'^^". These are 
common multiplicative factors, as are soft gluon re-summation effects, and cancel in the ratios taken. 

^ For other model-independent approaches to the determination of the Higgs couplings, see ifTOl [TTl l29l435l . 




(20) 



B. Global Fit To Higgs Signal Strengths 
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data [|36l - [39ll . and the recently reported Tevatron Higgs results BOl . The data we use is listed in 
the Appendix. We fit to the available Higgs signal-strength data, 

observed 



for the production of a Higgs that decays into the observed channels i = 1 ■ ■ ■ Nch- Here Nch 
denotes the number of channels, the label j in the cross section, aj^h, is due to the fact that 
some final states are summed over different Higgs production processes, labelled with j. The 
are the efficiency factors for the various production processes producing a final state j to pass 
experimental cuts. The reported best-fit value of a signal strength we denote by fii, and the we 
construct is defined as 



The covariance matrix has been taken to be diagonal with the square of the 1 a theory and 
experimental errors added in quadrature giving cTj. We necessarily neglect correlation coeffi- 
cients as these are not supplied. For the experimental errors we use ± symmetric la errors 
on the reported /ij, while for theory predictions and related errors we use the results of the 
LHC Higgs Cross Section Working Group [41]. The minimum (Xmm) is determined, and the 
68.2% (1 a), 95% (2 a), 99% (3 a) best fit regions are plotted as = Xmin + with the ap- 
propriate cumulative distribution function (CDF) defining the corresponding A^^. 

We first perform a one-parameter fit in terms of the free parameter Fg, that depends on the stop 
sector of NSUSY, and plot — xLin Fig- (left). This distribution is directly related to 
the fit in the broader space of the generic Wilson coefficients of the local operators contributing to 
/i — 7- 77 and gg — )■ h, as shown in Fig. [T] (right), although the match of the C.L. regions is only 
approximate due to the difference in number of degrees of freedom in the fit. It is not surprising 
that with the addition of this free parameter, the measure is improved as compared to the SM. 
However, it was by no means guaranteed that the stop line determined by the NSUSY relationship 
between the Wilson coefficients would pierce the best fit region away from the SM one. This 
accidental fact allows a ~ 2cr improvement of the fit. Although this is intriguing, we caution the 
reader that the interval of Wilson coefficients that intersect the 1 a best fit region is mapped into a 
very narrow range of stop mass parameters, corresponding to a tuned area of parameter space. In 
addition, best-fit regions in (c^, Cg) space more distant from the SM point at (0, 0) will generically 



Ei (Tj^h X Br(/i i)]sM 



(21) 




i=l 



(22) 
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Global Higgs Data 



7&8 TeV LHC & Tevatron data 
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FIG. 1: The one dimensional fit to Fg in the NSUSY scenario to global Higgs data (left), and the approx- 
imate projection of the relationship between the Wilson coefficients into the higher dimensional operator 
space (right).The green, yellow and gray regions correspond to the 1, 2, 3 o" allowed regions in the ID or 2D 
fit space (defined with the CDF appropriate to each case. This difference accounts for the mismatch in the 
Ax^'s that define the best-fit regions). Also shown as solid (brown) contours is the enhancement of the /x^-y 
signal strength and how such a condition projects into the best fit space. 

correspond to lighter states, as their impact scales as This will represent a further problem 

for this region. 

We can characterize the allowed relationship between the Wilson coefficients that intersect the 
(Icr) best fit region in a model-independent way, finding that current data is consistent with the 
following four ranges of the Wilson coefficient ratios, corresponding to the four different best-fit 
regions in that 2D space [28J. For > 0: 



-0.01 < Cg/c^ < 0.16 



0.27 < Cg/c^ < 2.5 



(23) 



and, for c^, < 0: 



- 0.1 < Cg/c^ < -0.065, 



-0.016 < Cgjc^ < 0.001 



(24) 



Note that these bounds are approximate in the following sense: for a ID fixed relationship between the Wilson 
coefficients, the allowed C.L. regions are slightly different if obtained with the ID CDF or for the 2D Wilson 
coefficient case. Again, this effect can be seen in the NSUSY case in Fig.|T| 
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FIG. 2: The projection into stop parameter space of the best-fit regions from a global fit to Higgs signal 
strength data. Here 5m = (m?^ — m? )^/^. Colour convention is the same as in previous plots. The inset 
zooms into the low mass best-fit region (ruled out by LHC monophoton searches). We have varied tan /3 
in the range (2, 20) and taken the overlap of the best fit spaces, slightly increasing the allowed parameter 
space. The three plots show the cases of no mixing {6^ = 0), intermediate mixing (9^ = it/ 12) and maximal 
mixing {6^ = 7r/4). The dashed line corresponds to the second minimum in the one parameter shown in 
the previous figure. 
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In the limit of a single field contributing to the Wilson coefficients, the Cg/c^ ratio is dictated by 
the quantum numbers of the field integrated out.^' Clearly, the study of the possible intersections 
of such lines with the best-fit regions in the space (cg, c^) for any model (including NSUSY) will 
become much more important with further refinements in the measurement of Higgs properties. 



We discuss some prospects for the improvement of these fits in Section VI 



For the NSUSY case, the light stop best-fit region occurs for Fg ~ 2; one can see how this 
space relates to the (c^, c^) plane in Fig.[T| and it corresponds to stops significantly lighter than in 
the region of small Fg. This space is already strongly constrained by monophoton searches, as we 
discuss further below. NSUSY hopes in light of current global Higgs data (when our assumptions 
are adopted) are based on the consistency of NSUSY in the (cg, c^) parameter space near the SM 
point (cg, c^) = (0, 0), for larger m^^ and small Fg. Translating the allowed fit space to the space 
of the stop parameters is very convenient to discuss the interplay with further constraints and direct 
stop discovery prospects. When we translate the results of the global fit to Higgs signal-strengths 
to the stop space, we find the best-fit regions shown in Fig. |2j The three plots show the cases 
of no mixing {Oi = 0), intermediate mixing (9^ = 7r/12) and maximal mixing (9^ = vr/4). We 
will show these canonical parameter choices throughout this paper when examining the global 
constraint picture. The point in the ID fit that has a local minimum with small Fg is mapped 
into the red dashed line in each of the stop space plots in Fig. [2] and, as expected, has larger m^^. 
As anticipated, the best-fit region at low stop-mass is extremely narrow: a sure indication of its 
fine-tuned nature. 



C. Br(S -). Xs 7) 

Another important challenge for the parameter space of NSUSY scenarios comes from non- 
SM contributions to magnetic moment operators. Although the contributions to these operators 
vanish P3l in the pure SUSY limit, NSUSY scenarios are far from this limit by construction. As a 
result, the reduction in the allowed parameter space due to constraints from Bi(B — )• Xs 7) can be 

See Ref. ll42l for a recent study that also emphasizes this point. 
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significant. Recall that the effective Lagrangian (neglecting light quark masses) is given by [|44l 



* 1=1 



Gf 



Vtb m, [C, SL a^'''hneF^, + C^-SLa^'-TahRg,Gl^ + ---]. (25) 



4\/2 7r2 

The SUSY contributions to the magnetic moment operators are well known [|45l and can be 
applied to the particular NSUSY scenario. The dominant contributions to the Wilson coefficients 
come from stop-chargino loops: 



2 



2 



— 1 3 si mi 

i=l I P h 



2 



- {U\tL){tR\ti) Fj^s 



2 






1 















(26) 



We have only retained the light Xi" with mass m ± ~ yuas consistent with NSUSY assumptions 



(the gaugino mass is M2 > yU ^ fnw) ■ The loop functions F^g are given in the Appendix. We 
vary the fi parameter in the range ~ 100 — 200 GeV. The lower limit of this range is set by LEP 
bounds on x"^ [|46l : the upper limit by naturalness considerations 161. In principle, there are also 
loop contributions from the light However, these contributions can be strongly suppressed in 
the minimal NSUSY limit we consider. For example, there is no source of breaking of the residual 
U(1)r symmetry, unless Majorana masses for the gluino are introduced. This symmetry plays a 
role in suppressing effects of large tan 13 or proton decay or flavour violating observables which 
require a flip in chirality, see Refs [|47ll48l . We will consider a minimal flavour violating scenario 
P9l453l when examining the NSUSY spectrum, this can follow from the argument above directly. 
In the case of Majorana gluino masses we assume MFV. 

We use the results of Ref. If54l for the constraints on the BSM Wilson coefficients Q = 
C^A/ _^ observable BR{B Xg 7) ij^>i.6 cev- The contribution of the BSM Wilson 

coefficients to this observable is given by 

BR(5 ^X,7)s^>i.6GeV = [(3.15 ±0.23) -8.0 AC7(/io) - 1.9AC8(/io)] x 10"^ . (27) 

Here we neglect (ACj)^ terms and have used (implicitly) the input values listed in Ref. [|54l . 
where, in particular, the scale fio = 160 GeV was chosen for the SM results. Another assumption 
used in Ref. [|54ll is that all other induced BSM operators have Wilson coefficients that satisfy 

^^|^~0(9U n>2. (28) 
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FIG. 3: Overlay of the best-fit regions of stop parameter space from global fit to Higgs signal strength data 
with the allowed space due to Br(i? — )• Xg 7) at the level of 1, 2, 3 o". The 3 a allowed region is overlaid 
with a vertical mesh, the 2 a has a horizontal and vertical mesh, while the 1 a allowed region is the green 
region with the further addition of the diagonal mesh. 
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and the NSUSY scenarios we are interested in will have to satisfy this condition when this con- 
straint is strictly applied. Comparing to the current world experimental average given in Ref. [|55]| 

BR{B ^ Xs 7)i?,>i.6GeV = [3.55 ± 0.24 ± 0.09] x 10-^ (29) 

we will use the 1(2) a bound -8.0 ACj - 1.9 ACg = 0.3 ± 0.34(0.69) to constrain the NSUSY 
parameter space. In Fig. [s] we show the interplay of the constraints from the global fit to Higgs 
data and constraints due to Br (5 — )• Xg 7) in minimal NSUSY. We see that there exists consistent 
parameter space that can pass both experimental tests, primarily at the level of ~ 2a in each case. 
Large mass splittings scenarios of the stop states when large mixing is present are significantly 
disfavoured. 



D. Electroweak Precision Data 



Measurements of mw and other EW precision observables also restrict the allowed parameter 
space of the NSUSY scenario. Recent measurements of mw at the Tevatron ||57ll2Hl are of partic- 
ular interest. The world average [|59l has been refined to {■mw)exp = 80.385 ± 0.015 GeV, with 
a significant reduction of the quoted error. As recently re-emphasized in Refs. [|60l [6TI . precise 
measurements of the value of mw constrain the allowed parameter space of a weak scale NSUSY 
spectra when nih is known. This occurs as the allowed custodial symmetry violation that could 
be present in the sfermion sector is bounded. Note that a global fit to EWPD produces a AT con- 
straint that has about twice the error of the constraint used here, which is directly determined from 
the shift in the W mass. 

In this Section, we add this further constraint in the study of the impact of NSUSY spectra on 
Higgs properties. We use the numerical approximation of the two-loop SM prediction of mw given 
by Ref. [f62ll and the method of Ref. [61] . The relevant SUSY correction to the SM prediction^^ 
of m is given in Refs. Il63] - l66l as 



(Amwr''''^ ^ , , AP"^^""- (30) 

12 We are aware of the recent Babai- result, BR{B 7)ij^>i.6GoV = [3.31 ± 0.16 ± 0.30 ± 0.10] x 10"*, 

Ref. Il56l . with a lower central value but larger error. As this result is not yet incorporated into a global average, we 
use the quoted HFAG average. 

'■^ Here sw is defined in the on-shell scheme. 
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FIG. 4: Constraints from precision measurements of mw, shown as 1,2,3 sigma shaded regions, over- 
laid on the allowed stop parameter space consistent with the global fit to Higgs signal strength data. Plot 
colour/display convention is the same as in previous figures. 



Neglecting terms proportional to small h mixing angles, the SUSY contribution is given by 
and the function Fq is defined as 

r 1 ^ *^ ^ 1 

Fo X, ?/ = X + y log -. 

x-y y 

We can trade mi in this expression for the stop masses and the stop mixing angle {6i) using 



m^^ ~ cos Oim^^ + sin Oim^^ — — m^y cos(2 (3). 



(31) 



(32) 



(33) 



This assumes small sbottom mixing and neglects perturbative corrections. We are also neglecting 
other SUSY contributions to /S.p^^^^ , namely the x^'^ contributions from the light states retained 
in the residual NSUSY spectrum. The masses of these states are set by the same scale (~ /i) in 
NSUSY, and the mass splitting of these doublets is small, leading to subdominant contributions to 
/\pSusY _ In order to define the SM prediction of mw we must specify a central value and error in 
nih- Given the current state of the data we simply take the Higgs mass range 125 ± 2 GeV for this 
constraint. The effect of this uncertainty in nih is very subdominant to the largest source of error, 
which is the uncertainty in (which leads to a ±5.4 MeV effect on m^y). 

With the choice of input parameters in Table |ll| we find {m]y)sM = 80.368 ± 0.006 GeV, 
which constrains NSUSY through the resulting bound (Amw)^^'^^ ^ 0.023 ± 0.016 GeV. This 
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translates into a constraint Ap^^^^ < (4.1 ± 2.8) x 10~^. This is in good agreement with the 
result of Ref. ll6T1l . which uses the same method, up to small differences in the input parameters. 
Note that for Oi> n /A the lightest stop is dominantly In compared to the interval we discuss, 
< 6i < n / A,'m which it is dominantly II- However, in the former case one still obtains a similar 
constraint space for the /S.mw constraint, with the shift in the SM prediction and the best fit value 
selecting for degenerate stop states. 



1. The Funnel Region 

We show the overlap of the global Higgs fit constraints and the constraints due to mw in Fig.|4} 
Note the good degree of consistency between both constraints. This consistency can be traced to 
the following. The low mass "funnel region" in stop parameter space arises from minimizing the 
contribution to Fg. The condition — )■ translates into the relationship 

2mt (m'j + m? )^/^ 

sin(2 6^) ^ (34) 

(mj — mi Y 

Using this constraint in the maximal mixing case, and substituting into the Apg^"^^ expression, 
one finds (expanding in mt/mi^ < 1) 



_^ m 



2 m? 

12 



(35) 



This corresponds to a suppression of Apg ^"^^ to the right order of magnitude to match the condi- 
tion Ap^^^^ ~ 10^"^. This basic correlation also occurs outside of the maximal mixing case we 
have considered for simplicity, although the above equation is modified by 0(1) numerical coef- 
ficients. Conversely, when the stop parameters are not related in this manner, Apg^'^^ receives a 
contribution dominated by the largest mass scale 

/^o 8v^7r2^ *^ '''' 

and is significantly larger, ruling out the associated parameter space. In this sense, this region is 
clearly associated with some cancelations due to parameter tuning. 



E. Collider Bounds 



In this Section we describe some of the current collider bounds on the NSUSY spectrum and 
how they pertain to the stop parameter space of interest found in previous sections. The most stud- 
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ied and stringent bounds on NSUSY come from missing energy signatures. However, including 
R-parity or not in a weak scale NSUSY spectrum is not dictated directly by naturalness. Due to 
this, we will mostly restrict our attention to more robust collider constraints. We note however 
that with the absence (to date) of missing energy signatures, several searches for R-parity violat- 
ing phenomenology (without significant missing energy) [18] in the multijet and multilepton final 
states are now ongoing. These studies could further constraint the interesting parameter space that 
we isolate in this study. 

The main features of NSUSY relate to the Higgsino and Stop sectors. As we have discussed 
before, if Higgsinos get their mass solely from the ji term, we expect a very small mass splitting 
between the x^? of the order of /x/Mf 2)- Here Mi 2 are the gaugino masses that are taken 
to be large Mi 2 ^ in minimal NSUSY. In the deep-Higgsino region, the — )■ x° + decay 
occurs with no hard-pr objects to tag on. Then, only monojet/monophoton searches at LEP were 
sensitive to such decays, leading to a bound m^^i > 92 GeV [|46ll . 

Searches for stops are limited by the complexity of the final state and the similarity with the 
SM tt background, especially in a scenario where the chargino and neutralino are heavy. The 
recent direct stop combination searches from ATLAS [j67l do not have sensitivity on the stop 
mass, except in a very small region around nii ~ 300 GeV, and only for a limited range of 
One can obtain more information on the stop sector by looking at searches of a single photon 
plus missing energy [|72l 1731 . Although those searches were performed in the context of large 
extra-dimensions and dark matter effective theories, a straightforward re-interpretation in terms 
of SUSY can be done 11741 . within the assumption that the stop and the neutralino are separated 
by < 30 GeV. This is indeed the case in the low-mass region preferred by the Higgs fit, namely 
nil < 120 GeV. Since the Higgsino is heavier than 90 GeV, the splitting between x and t is right 
at the best sensitivity point. The mass bound obtained in Ref. [|74l is 150 GeV (95% C.L.). This 
strongly constrains/excludes the low mass stop parameter space most preferred in the global Higgs 
fits. Monojet searches could also be used in the low m^^ region of parameter space 117511 , and one 
could also use top precision measurements (spin correlations) to rule out this area of parameter 
space li69l . 

Finally, we note that gluino-assisted stop production searches rely on a light gluino around a 



Note that there are some proposals to improve the stop searches Il68ll7n . using different strategies, from boosted 
tops to fitting the missing energy distribution. 
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TeV, which is in some tension with the combination of multijets+leptons+missing energy [|76ll77]| . 
which leads to a bound of 1050 GeV. This gluino bound is rather model independent, as inter- 
mediate x^, X2 to leptons or off-shell stops do not change the bound considerably fTSl |79|.^^ 
Nevertheless, a gluino heavy enough to disable this search is allowed in a minimal NSUSY sce- 
nario and the stop limits from these gluino-assisted studies do not (as yet) directly constrain the 
parameter space of interest. 

V. COMBINING CONSTRAINTS 

It is of interest to address how an NSUSY scenario globally fits the increasingly rich dataset. 
In past sections, we have shown the interplay of various precision constraints; in this Section we 
perform a global fit to Higgs signal-strength data, constraints due to precision measurements of 
m\Y and Br(i? Xg^). We fix tan /3 = 10 and show in Fig. [s] such global fit results for the three 
mixing cases we have considered in this paper, to illustrate the allowed global fit space. 

The results in Fig. |5] show a good fit to the data. We are fitting to fifty observations: forty-eight 
Higgs signal-strength measurements, as well as /S.mw and Br (5 — > 7). For the three mixing 
cases 9i = 0, 7r/12, 7r/4 and /i = 100, 200 GeV we find, fitting to the two stop parameters m^^, 5m, 
that Xmin/ (50 — ?)) ^ 1. The interplay of the constraints that allows a good fit is of interest. The 
best-fit region in the Higgs signal- strength fit at very low masses, m^-^ ~ 100 GeV, is essentially 
ruled out due to mono-photon searches and constraints from Amvi/ and Br (5 — )• Xs7). When 
combining constraints, one finds a best fit to the Higgs signal- strength data with a larger m^-^ (with 
less mass splitting for larger 61) in a manner that also improves agreement with the small deviations 
from the SM predictions in Amvi/ for large regions of parameter space. However, the significance 
of this observation is currently marginal: the SM gives a comparable fit to these observables with 

xLjm ^ 1-0. 

A. Fine-tuning and the Higgs Mass 

It is also interesting to examine if the Higgs mass could be consistent with its observed value 
in the allowed parameter space of the fit, without large fine-tuning. We have included the large 

Although this search is optimized for Majorana gluinos, the Hmit would be even more stringent for Dirac 
gluinos Ho). 



21 




200 400 600 800 1000 200 400 600 800 1000 200 400 600 800 1000 




200 300 400 500 600 700 800 200 300 400 500 600 700 800 200 300 400 500 600 700 800 



^t, 

FIG. 5: Best-fit regions in stop parameter space from the global x^- Colour conventions the same as in 
previous figures. We show the three mixing cases 6^ = 0, 7r/12, 7r/4 for two values of fi that set the 
Higgsino mass scale in these minimal scenarios, /x = 100, 200 GeV. The tan /? dependence of the result 
is primarily driven by the Br(^ — )• Xs 7) constraint. Larger values of tan /3 select for a more degenerate 
stop spectrum when 6^ 7^ 0. In determining Xmm' '^^ have assumed a prior 150 GeV < , due to 
monophoton constraints. In each figure the best fit point is marked with a dot. 

radiative corrections to the Higgs mass expected after SUSY breaking [81] by using the public 
code FeynHiggs [[82l [83l (for alternatives, see [[841 '). which includes all one-loop corrections and 
the dominant two-loop effects [[85 fl . This gives a precise enough determination of rrih (with an error 
of a few GeV), provided the hierarchy in the stop masses is not so large that further re-summation 
of logarithms is needed [[861 . Precision is needed here because the soft masses required to give a 
large enough Higgs mass are exponentially sensitive to ruh. 

We show in Fig. |6]the band in stop parameter space that gives nih = 125 ± 2 GeV, overlaid 
on the best fit regions. The Higgs mass has a non-negligible two-loop sensitivity to the value of 
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FIG. 6: The overlay of the Higgs mass condition mh = 125 it 2 GeV (in red/darker shaded region) and 
the best fit space (colour convention the same as in previous figures). Also shown are labeled contours of 
constant fine-tuning (soUd black) for a UV cutoff A = 10 TeV. The same contours are shown (dashed grey) 
for A = 50 TeV. In the no-mixing case, the Higgs mass constraint cannot be accommodated in minimal 
scenarios: no Higgs mass band is present in these plots. 

the gluino mass, which we have chosen to be 1.1 TeV in this figure.'^ For non-zero stop mixing, 
larger values of the gluino mass tend to lower rrih and therefore require larger values of the stop 
masses [|82l . We see in Fig.|6} left plot, that the best- fit region with sub-TeV stop masses and zero 
mixing is not consistent with the Higgs mass condition, as was to be expected. Larger stop masses 
with large mixing can however produce the observed Higgs mass value (Fig. |6} middle and right 
plot) in exceptional regions of parameter space (or when a worse fit is considered), although the 
degree of fine-tuning in these parameter regions is a concern. 

It is well known that physics beyond the minimal supersymmetric version of the SM can easily 



We take the other sfermion states to be ^ 2 TeV in this analysis, however the Higgs mass band is insensitive to the 
particular mass values of these states for tan /3 ^ 10 and sfermion masses in the TcV range. 
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FIG. 7: Effect of an upward shift on the Higgs quartic, 6X on the allowed Higgs mass band rrih = 125 it 2 
GeV. We have chosen the stop parameters as for one of the cases of the previous figure, see inset of lower 
right plot. The corresponding Higgs mass band is the small red region on the right. For 6X = 0.005 (0.01) 
this region extends to lower stop mass values and is given by the faded yellow (orange) band bounded by a 
dashed (dotted) line, as indicated. 

give contributions to the Higgs quartic coupling which controls the Higgs mass. This can happen 
through new sectors coupled to the Higgs in such a way that they add at tree-level new supersym- 
metric F-term or D term (or even susy-breaking) contributions to the Higgs quartic even in the 
Higgs decoupling limit (see ||87] - [89|| for an incomplete list of studies on this). A very moderate 
tree-level upward shift of this quartic, which can be interpreted as a threshold correction at the 
SUSY scale, would have a direct and very important impact on increasing the predicted Higgs 
mass and reducing the associated fine-tuning. We have illustrated this in Fig. |7| which assumes 
an additional contribution 5X = 0.005,0.01 to the Higgs quartic coupling A = m^/(2t>^). The 
new physics required for this shift can, on the other hand, have a negligible impact on the cou- 
pling of the Higgs to gluons or photons. Therefore, in such scenarios one is allowed to soften the 
link between nih and the stop sector, while the analysis of the impact of Higgs search results on 
constraining stop parameters we have performed would still apply. 

It is difficult to define a uniquely compelling fine-tuning measure in an effective theory, or argue 
what degree of fine-tuning is clearly unacceptable. Fine-tuning considerations are necessarily 
dependent on the UV physics of the EFT and so one can only make a rough estimate in a low- 
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energy effective theory. We could follow the definition of fine-tuning measure of [19011911 . which 
quantifies the tuning by measuring the sensitivity of the electroweak scale (as given by e.g. the Z 
mass) with respect to changes in the fundamental UV parameters. However, as we want to keep 
an open mind about what detailed UV physics completes the NSUSY scenario, we will content 
ourselves with a different estimate of tuning that simply compares the value of the Z mass with 
known loop contributions to it (in our case those coming from stop corrections), requiring that the 
latter are not much bigger than m^. That is, defining 



A; 



(37) 



the associated fine-tuning will be 1 part in A^. Restricting ourselves to the moderate value of 
tan /3 ~ 10 and assuming that stops loops are the dominant contribution to this fine-tuning mea- 
sure, one finds (see e.g. [|92ll ) 



m? + m? 

Il 12 , 



3 



;7r2t;2 



2A2 

log I 2 ^ 2 I • (38) 



The scale A is associated with new states required to cut off the logarithmic divergence in the 
effective theory, and is associated with the messenger scale that transmits SUSY breaking from 
a hidden sector. The choice of this scale is UV dependent, offering some further caution on the 
interpretation of the results. For numerical purposes we consider this scale to be A ~ 10, 20 TeV 
consistent with recent choices in the literature [,6J. The fine-tuning contours are overlaid in Fig.|6] 
in this case (for related work, see Ref. [|93l ). 

It in interesting to note the exceptional region of parameter space in the global that allows 
the Higgs mass to be raised with relatively light stops and smaller fine-tuning, see Figs. |6] and |7] 
for = 200 GeV (light stop mass insets). This is another example of the "Funnel region" of stop 
parameter space. The Br (5 — )• Xg 7) constraint is sensitive to this region of parameter space, and 
we note that although currently allowed with the HFAG average we have used, incorporating the 
recent Babar update into the a global average of Br(i? — )■ Xg 7) will further resolve this interesting 
region of parameter space. We also note that, if this parameter space continues to offer a good fit 
to the data, such stop masses are very difficult to probe directly in collider searches. A continued 
deviation consistent with such stop parameters in the global Higgs fits could be the leading exper- 
imental indication of the presence of stop states in this exceptional region of parameter space. We 
discuss the prospects for such future fits in the next Section. 
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FIG. 8: Projected best-fit areas when 8 TeV signal-strength values hypothetically report by the end of 2012 
have their errors scale down to ai/2 while they keep their current 8 TeV best-fit central values (left) and 
when the best-fit central values converge to the SM (right). For comparison, we also show in faded gray 
(solid, dashed, dotted) lines the current (1, 2, 3) cr best-fit contours. 

VI. PROJECTIONS FOR NSUSY AT THE END OF THE 8 TEV RUN 

A. Fit prospects. 

The integrated luminosity of experiments at LHC is increasing at a rapid pace, and the power 
and precision of global Higgs fits of this form have scaled remarkably to date in the 7, 8 TeV 
runs. If this scaling continues (i.e. if systematics and/or correlations do not become the dominant 
source of error by the end of this year) and each experiment gathers the projected ~ 30 fb^^ before 
shutdown, we can study the NSUSY prospects for fits of this form by the end of the 8 TeV run. 
We reduce each of the reported 8 TeV signal strength errors by a factor 2, down to ai/2, and 
assume two different scenarios for the final 8 TeV fii central values: 1) they retain their current 
central values, or alternatively, 2) they move to /ij = 1. As the integrated luminosity is to increase 
by more than a factor of four while more channels are also expected to be added to the global 
Higgs dataset, this estimate is conservative. We show the prospective best-fit regions in the two 
dimensional Wilson coefficient space of the hgg and h'yy operators in these two hypothetical cases 
in Fig. [8] It is clear that a much more constrained stop space due to these Higgs fits, and possibly 
a 1 — 2a exclusion of NSUSY, will be feasible by the end of the 8 TeV run. Due to decoupling, 
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if the SM emerges as a better fit in the dataset, with /tf ^'^^ 1, NSUSY will be relatively less 
constrained by fits of this form, however, fine-tuning arguments could be used against it. 

B. Prospects for bounding the best-fit regions: collider and Higgs search exclusions 

In this Section we compare the impact of the global fit on the stop parameter space with the 
searches based on missing energy. The global fit we have performed with NSUSY draws interest- 
ing conclusions: light and nearly degenerate stops are preferred in the range of 400-500 GeV and 
there is an exceptional region of parameter space with even lighter stops. These regions accommo- 
date the global Higgs properties and the mw and Br (5 — )• Xs 7) constraints. On the other hand, 
if the value of the Higgs mass is due solely to the stop sector, the allowed region shifts. Somewhat 
heavier stops, with a larger separation (> 200 GeV) seem to be indicated in this case, mostly out- 
side of the exceptional region of parameter space. The same fit to Higgs data does rule out a small 
portion of the stop parameter space, comparable to the monophoton exclusion. Note however, that 
those two exclusions, from monophoton searches and from Higgs data, have different degrees of 
SUSY model dependence. The monophoton exclusion is based on missing energy signatures due 
to assumed R-parity conservation, whereas the Higgs data exclusion does not assume it, or any 
particular stop decay chain. 

In the experimental searches, stops are assumed to decay to In the case of NSUSY, those 
electroweak states are degenerate higgsinos, leading in both cases to missing energy in association 
with a b-jet or a top, 

i — > 6 x"*" ~^ + soft objects 

i^tx'. (39) 
The first decay chain corresponds to the direct sbottom search [|60l |94] - |97ll . The second decay 



contains a top (on or off-shell), considered by the direct stop searches by ATLAS (See Sec. IV E I. 

The current stop searches are very weak in the region of m^n ~ > 90 GeV. The current 
ATLAS exclusion region is a small triangle below /i =160 GeV, and between, roughly, 300 and 
450 GeV. To explore NSUSY in this best-fit region, one would need to push the stop searches 
not towards larger values of m^, but of m^o. The main issue to reach higher values of fi is the 
similarity with tt when either for nif rrit + n, or large /x. In both cases, the missing energy 
distribution loses its ability to suppress the top background. The stop has little phase space to 
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FIG. 9: Left figure: The missing energy distribution for different choices of neutralino mass, with the stop 
mass fixed to 450 GeV. The red-solid, black-dotted and blue-dashed Unes correspond to m^=100, 200 and 
250 GeV. Notice the decrease of missing energy as the neutralino mass increases. Right figure: Restriction 
in the stop space due the limits on the mass splitting of 6/,, t^, with a lower bound on nij^ > 200, 400, 600 
and 800 GeV. The vertical lines correspond to interpreting direct sbottom searches as a t — )• bx^ decay, 
leading to the same mass limit. 

produce boosted x"^'^ unless it itself comes with some boost. For ~ 300 GeV, this would need 
some requirement on radiation jets. We illustrate the reduction of missing energy as one increases 
the value of fi in Fig.|9j Similarly, direct detection of charginos is difficult, due to the degeneracy 
with the neutralino and small cross sections. 

Direct sbottom production searches can be translated into stop production searches when the 
stop decays to b and x^, and this search is more encouraging. The ATLAS collaboration performed 
a search for sbottom squarks, resulting in a 95% C.L. upper limit m^^ > 390 GeV for neutralino 
masses below 60 GeV with C = 2.05fb^^ [|98l , but would probably be extended to larger values of 
m^o with the 2012 dataset, being able to cover the case m^o,± < 200 GeV. Note that direct sbottom 
searches assume BR(6 — )• bx^)=\ but can be re-scaled for lower branching ratios. Moreover, 
pushing the sbottom limits, within NSUSY, is also an indirect probe of the stop sector due to the 



custodial relations shown in Sec. IV D This interplay between direct and indirect limits on stops 



from direct sbottom production searches is illustrated in Fig. 1 1 , where a maximal mixing case is 
shown. Interestingly, if the lightest stop is L-dominated, the combination of sbottom and custodial 
violation limits could be more constraining than the re-interpretation of direct searches in terms 
of stop decays into charginos. As we do not know the chiral admixture of the stops and sbottoms. 



28 



using both strategies would give the best sensitivity. Since sbottom direct detection searches do not 
suffer many of the problems that affect the advance in collider stop searches, this further supports 
the observation that direct sbottom searches may be the best direct access to stop parameter space 
in the near term [[60l . 



7&8 TeV LHC data & Tevatron 




FIG. 10: The region of the higher dimensional operator space excluded at 95% C.L. from current Higgs 
searches is shown by the black shaded area overlaid on the best-fit regions. The NSUSY relation between 
the Wilson coefficients is shown with a red line. 



Higgs search data also provides a powerful insight into stop parameter space. In particular they 
offer experimental reach into the large jj, region that is such a challenge for collider searches. In 



Fig. 10 we show the region of the higher dimensional operator space excluded at 95% C.L. from 
current Higgs searches. Translating this exclusion into the stop parameter space does not lead in 
general to direct lower mass bound that is independent of 5m (when 6i ^ 0) due to the existence 
of the "funnel region" where a cancelation of the stop contributions to Fg can occur. However, by 
jointly imposing the condition that the Br(i? — t- Xg 7) constraint is within its 2 a allowed region, 
we can define more stringent current exclusion regions for approximately mass degenerate stops. 
By studying how the Higgs data will scale with more luminosity, we can also project expected 



exclusions in the stop space by the end of the 8 TeV LHC run. We find the results shown in Fig 1 1 



In the case of these joint exclusions, there is in general a tail, or spike, of allowed parameter 
space in the Br(i? — )■ Xg 7) constraint where a significant degree of cancelation is occurring in the 
NSUSY contributions to Br(i? — )■ Xg'^). This reduces the scaling of this bound in general. This 
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feature can be seen in Fig. |3] (bottom, right) and the viability of this cancelation region is sensitive 
to further revisions of the global experimental average of Bt{B — )■ Xs 7) including the recent 
Babar measurement. If one does not rely on such cancelations to diminish the interplay of the 
Br(i? — 7- Xg'j) and Higgs search data constraints (outside the funnel region) then the exclusion 
bounds scale much more quickly, as is also shown in Fig 1 1 

Note that these exclusions are in the region of /i ~ 100 — 200 GeV in NSUSY and do not 
require a missing energy tag. This makes indirect studies of stop exclusions from Higgs search 
data broadly applicable to many SUSY models and these bounds have significant reach into the 
large /i region. Conversely, collider based searches will be severely challenged to reach the large 
/i space in direct stop searches due to the fact that the signal region is increasingly similar to the 
t i SM backgrounds with small missing energy. 

VII. CONCLUSIONS. 

In this paper we have examined the light that Higgs search data sheds on the stop sector of 
Natural SUSY, where stops are expected to lead to the largest effects on Higgs properties. We 
have interpreted a generic fit to the latest Higgs data in models which allow for additional BSM 
contributions to the hgg and h'j'j couplings in terms of stop parameters. In performing this careful 
study of the impact of light stop states on Higgs properties, we have included QCD matching 
corrections and consistently treated the allowed best-fit regions with a one parameter CDF. Further, 
we have examined in detail the impact of related precision flavour and EW constraints. 

Interestingly, we find that a combined global fit to Higgs properties, Bi(B — > Xs 7) and Amw, 
show a preference for ~ 400 — 500 GeV degenerate stops, that is largely independent of the stop 
mixing angle. Such a spectrum of stop states cannot explain directly the Higgs mass value of 
~ 125 GeV but they do offer a good fit to the global dataset we have studied. The quality of 
fit is comparable to the goodness of fit in the SM. Such a stop spectrum can also improve the 
agreement of the observed and predicted mw, ameliorating a very slight tension in the SM with 
this observable that has arisen with recent precise measurements at the Tevatron. However, we 
emphasize that what is more interesting at this time is the degree of consistency that the small shift 
in this precision observable has with the allowed low mass region in the Higgs global fits. 

Further, we have also uncovered an exceptional region of parameter space with a light stop 
eigenstate m^^ ~ 220 — 260 GeV and larger mixing, with a second eigenstate that is significantly 
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FIG. 11: The current (future) exclusion regions from current (and prospective) Higgs search data that are 
directly combined with the condition that the Bt{B — )• Xgj) constraint is within its 2 a allowed region. 
Shown are the summary 95% C.L. exclusion results from the dedicated ATLAS collider searches for direct 
stop production in the m^^ , x? plane as solid regions (generally below the line x? = 90 GeV). These results 
are reported in Ref. |f67l| and were combined into the presented form at ICHEP2012. Added to these figures 
in the x? — 90 — 200 GeV space are current exclusions based on Higgs data and the Bt{B — )• Xs 7) 
constraint as a green region with diagonal lines; future exclusions from Higgs search data that evolves to 
more SM-like values (blue region with horizontal lines, labelled "SM-like") or when the global Higgs data 
has evolved to have the current 8 TeV central values while its errors have been reduced by a factor of 2 (red 
region with vertical lines, labelled "77 enhanced"). Top (bottom) figure is the no (maximal) mixing case. 
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heavier. In this exceptional parameter space region, the presence of such stop states can offer a 
good global fit to the dataset we have studied of quality comparable or better than the SM and 
can raise the Higgs mass to its observed value in NSUSY with fine-tuning ~ 1 — 5% (according 

to the measure we have adopted). Such stop states can also ameliorate the slight (statistically 
insignificant) tension in the SM prediction and measurement of mw- Future careful studies of this 
form will be essential in ruling out, or confirming a discovery of Natural SUSY at the LHC. 
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Appendix 
1. SM Inputs Used 
2. Efficiency corrections due to Higher Dimensional Operators 

We seek to draw as precise a conclusion as possible about the Higgs signal strength fit to higher 
dimensional operators with Wilson coefficients Cg, Cj in this paper. The effect of these higher 
dimensional operators on the signal strength parameters can have two forms. Directly the o-gg^h 
and Th^^^ rates can be effected, as characterized by our rescaling the data with the effects of Cg, Cj 
consistently. This effect we fit to directly using our global fit procedure. 

There is also a further effect that higher dimensional operators can have on the event yields 
that lead to the /Xj. The higher dimensional operators can alter the shape of the final kinematic 
distributions (that sum over more than one production mechanism), leading to a further correction. 
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Input 


Value 




125 ± 2 GeV 


mt 


173.2 ± 0.9 GeV El 


m;,(lS) 


4.65 ±0.03 GeVniOOlllOll 


rric 


1.275 ±0.025 GeV niOll 


rrir 


1776.820.16 MeV llOl| 


mz 


91.1876 ±0.0021 GeV lIlOll 


as(Mz) 


0.1184 ±0.0007 GeV |102| 




(275.7 ± 1) X lO"'' [103 ] 




(314.97686) x 10"^ |l04l 




0.2233 ± 0.0002 lIlOU 




1.1663787 ±0.0000006 x 10"^ GeV-2||101| 



TABLE II: Input values used in determining constraints on the NSUSY parameters in the indirect tests. 
When we use the value as{mt) we determine the value from NLO running in QCD finding the central value 
as{mt) = 0.1080. 

as a function of the Wilson coefficient, on the signal strength parameter. The Cg, do not affect 
the shape of the kinematic distributions due to of an individual production process, such as gg 
/i — > 7 7, however, they do affect the relative proportions of gg versus Higgstrahlung and YBF 
initiated Higgs production. This later effect is an "efficiency correction" on the change in the 
number of expected events expected due to a different effective efficiency for passing the cuts 
of the experimental analyses when higher dimensional operators are present. This effect is not 
captured in our direct fit to Cg, and should be quantified for precise conclusions. 

We have examined the effect of these efficiency corrections through simulating the effect of 
those operators in the 7 TeV run data in the — )■ /i — t- 77 (Oj, Ij, 2j). We added the shift due to 
the operators in the hgg and h'yy effective vertices using Feynrules [I105II into an UFO model [I106II 
of MadGraphS H 10711 . We generated samples at 7 TeV with a parton level cut on \r]j\ < 3.5 and 
Pt,7 > 20 GeV. The parton level events generated by MadGraphS are passed to Pythia II108II to 
simulate the effects of parton showering, and then to Delphes Ill09ll for a fast detector simulation. 
We use the generic LHC parameters for Delphes, and reconstruct jets with the anti-fcr algorithm 
using 0.5 for jet cone radius. We cannot simulate all the characteristics of photons in the different 
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bins, especially the converted/unconverted nature of the photon. Instead, we check the effect of 



the new physics on the basic selection cuts described in Table III Those are on top of the detector 
effects simulated by Delphes, including isolation and energy-momentum smearing. 



ATLAS cuts 


CMS cuts 


^71,2! <2.37 


IVfiJ < 2-5 


> 40 GeV, > 25 (30) GeV 




m^^ = 125(126) ± 3 GeV 


m^^ = 125 ± 3 GeV 



TABLE III: Cuts implemented in the Monte Carlo simulation after Delphes. Numbers in parenthesis corre- 
spond to the 8 TeV run. 



Those cuts are then applied to the signal, containing the SM and the new operators. In Fig. 12 
we show the effect of the operators, c^,Cg. The effect on the cross section is sizeable, whereas the 
effect on the efficiencies due to basic px, t] and Am^^ cuts is moderate (less that 1%) except when 
the cross section drops due to a cancellation between the SM and new physics contributions. 
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FIG. 12: Left: The cross section as a function of Cj, Cg using the fixed relationship between the Wilson 



coefficients. Right: Change in efficiencies to the cuts defined in Table. Ill Both figures correspond to t/s 
7 TeV. 



The effect on efficiencies for the relation between and Cg given in Eq. 20 is very similar to 



Fig. 12 namely of the order ^ 1 %. Further the corrections stay at this level when both operators 
are present or when the effects at 8 TeV are simulated as we have explicitly verified. Due to 
this small correction we neglect effects due to efficiency corrections due to higher dimensional 
operators in this paper. 
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3. Br(S — >• Xs 7) Loop Functions 

The loop functions in Br (5 Xs 7) are given by Il45l [541 [TTOl 

, x(7-5x-8x^) x2(3a;-2), _ 
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All Data (48 channels) 



WW^Uw : 0.32±0.83 
bb : 1.97±0.71 
yy : 3.62±2.76 
TT : 0.45±1.8 
bb : 0.46±2.18 
ZZ^41 : 1.11 ±0.92 
WW^llvr : 0.45±0.62 
7)'uciPTi:''-54±1.46 
rruchPT,:0.21±1.98 
TTuRiPTi : 2.52±1.69 
r/uRhPT, : 10.42±3.68 
TYccivrt '■ 6.08±2.63 
rrcchPTi : -4.36±1.78 
TTcRiPTi : 2.74±2.03 
rrcRiiPTi : -l-63±2.87 
77ct : 0.35±3.61 
TTfi : 2.71±1.9 



ATLAS 7 TeV 



ZZ^41 : 1.03±0.64 
WW-.IIW : 1.86±0.65 
TTuciyr.: 1-W±1.3 
ryuchFn:0-89±l-74 
7>'m,pT,:0.95±1.39 
ryuMPT, : 1.88±1.84 
r7cciFTi:3.48±1.95 
>TcchPTi : 3.57±2.59 
rrcRiPTi : 0.47±1.74 
yycRhFri:0-4±2.1 
jY^j : 5.55±3.3 
yj'j,j2^58± 1.7^3 ^ 
rr:0.62±1.13 
bb : 0.59±1.17 
ZZ^41 : 0.64±0.57 
WW->llw : 0.38±0.56 
77„ : 3.16±1.81 
77, : 0.66±0.95 

772 : 0.74±1.16 

773 : 1.53±1.6 
77:: : 4.21±2.05 



ATLAS 8 TeV 



TT : -0.72±0.97 
bb : 0.41±0.94 
ZZ^41 : 0.79±0.56 
WW^Ilyv : 0.98±0.71 
77„ : 1.45±1.25 
)'7, : 1.5±1.04 
77, : 0.95±1.16 
77, : 3.76±1.77 
rr„i„o>= : -0.61±2.03 
yri„iEh, : 1-32±1.58 
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FIG. 13: Summary of Higgs best-fit signal- strengths /xj used in our global fit. For more details see Ref. f26l- 
l28l . Shown are the reported fii (or the reconstructed 8 TeV value if not directly reported II28II ). For CMS 
and the Tevatron we use values at m/j = 125 GeV, while for ATLAS we use m/j = 126.5 GeV. This choice 
is partially due to the limited experimental information currently supplied. For discussions on possible bias 
from this choice see Ref. ll28l llllll . Also shown is the combined jl and error as a vertical green band and 
the SM expected signal strength as a black vertical line at /i = 1. 
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